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We study the phase diagram of spatially inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov super- 
conducting state using the Ginzburg-Landau free energy, derived from the microscopic hamiltonian 
of the system, and notice that it has a very clear Lifshitz tricritical point. We find the specific heat 
jumps abruptly near the first-order line in the emergent phase diagram which is very similar to the 
recent experimental observation in layered organic superconductors. Comparison with experimental 
data allows us to obtain quantitative relations between the parameters of phenomenological free 
energy. The region of the phase diagram where the specific heat jumps can be probed by doing a 
dynamical analysis of the free energy. 



PACS numbers: 74.25.Dw, 74.20.De, 65.40. Ba 

Almost 50 years ago Fulde and Ferrell[T] and Larkin 
and Ovchhinikovf^ made a startling prediction about 
the possibility of a spatially inhomogeneous supercon- 
ducting state at high magnetic field and low tempera- 
tures. The new-proposed state was unique as it had 
a spatially-modulated order parameter, in contrast to 
the spatially-homogeneous order parameter of a standard 
Bardeen-Cooper-Schrieffer superconducting state. This 
Fulde- Ferrell-Larkin-Ovchinnikov(FFLO) superconduct- 
ing state, other than being theoretically interesting, has 
major technological importance due to its high supercon- 
ducting current densities and thus remains a very active 
field in both theoretical and experimental research till 
date. 

The FFLO state, from a theoretical perspective, is 
quite intriguing as it retains its superconducting prop- 
erty overcoming the orbital and Pauli-paramagnetic pair- 
breaking effects, even at very high magnetic fields. This 
state has been vigorously studied for different physi- 
cal systems [SHS] which vary from heavy- fermionic su- 
perconductors to dense quark matter. The experimen- 
tal search for the FFLO superconductor is all the more 
interesting as, due to the stringent experimental condi- 
tions needed, the conclusive evidence of FFLO state is 
still lacking. The FFLO state will only occur if the su- 
perconductor is in the clean limit, i.e. the mean free 
path is much larger than the coherence length and the 
Maki parameter, a measure of the relative strength of or- 
bital and Zeeman coupling to the external magnetic field, 
a = ^/2Horb/ Hp is greater than 1.8. The most promis- 
ing candidates obeying these conditions are the heavy- 
fermionic superconductor CeCoIua [HHS] and layered, or 
quasi two-dimensional(2D), organic superconductors. 

The layered organic superconductors have become the 
focus of late as they are, in most cases, clean su- 
perconductors and the orbital pair breaking effect is 
mostly suppressed when a strong magnetic field is ap- 
plied parallel to the conducting layers. In an ex- 



periment done on the layered organic superconduc- 
tor /t-(BEDT-TTF)2Cu(NCS)2, m which BEDT-TTF is 
bisethylenedithio-tetrathiafulvalene, Lortz et al[14] first 
provided a true thermodynamic evidence for an increase 
of the upper critical field and showed the existence of a 
narrow intermediate superconducting region, which they 
claimed to be a FFLO state. This claim had been further 
supported by recent experiments on the same material 
while looking for magnetic torque properties [T3] and in 
nuclear magnetic resonance experiments [TB]. One of the 
most prominent features marking the onset of the FFLO 
state is an anomalous jump in the specific heat|14l fT7] . 
Though some very recent papers by Buzdin et al[18l [19] 
report a timely study on anisotropic effects in layered su- 
perconductors, the specific heat anomaly remains largely 
unexplored. 

In this paper we present a detailed study of the 
FFLO phase diagram with the help of Ginzburg- 
Landau(GL) free energy, calculated from the microscopic 
hamiltonian|20|. We calculate the specific heat to show 
that, indeed, the onset of the FFLO phase is accompanied 
by a jump in the specific heat across a first order line. 
This, we show, is a direct consequence of the fact that 
the FFLO phase diagram invariably exhibits a Lifshitz 
tricritical point(LTP) which is an outcome of keeping 
the relevant, in renormalization group terminology [21], 
higher order gradient terms in GL free energy. The region 
in the phase diagram where the discontinuity in specific 
heat occurs can also be experimentally ascertained from 
a study of the temporal behaviour of dynamical struc- 
ture factor. We have shown this by suddenly quenching 
the system from the normal to the FFLO state across 
the LTP and calculating the structure factor using the 
time-dependent Ginzburg-Landau(TDGL) formalism. 

We begin by considering the action of a two- 
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component, spin-imbalanced Fermi gas 

+ g J dT(fx{ij'^ijY. (1) 

■0 = i'^UTi'd)^ is a column matrix containing the an- 
nihilation operators of spin-up and spin-down fermions, 
respectively, and r is the imaginary time. As we are 
looking for the finite-temperature properties of this sys- 
tem, we need to use finite-temperature(or imaginary- 
time) Matsubara formalism, fj, + Sfj, and fx — Sjj, are the 
chemical potentials of the up- and down-spin fermions, 
respectively. Clearly, 5pL is the measure of imbalance. 
In a superconducting material this spin polarization re- 



sults from the Zeeman coupling of the electron spin to a 
magnetic field, as discussed earlier. In ultra-cold, imbal- 
anced Fermi gases mismatch in spin polarization causes 
this chemical potential difference, "g" is the attractive 
electron-phonon coupling constant for layered supercon- 
ductors and for an imbalanced ultra-cold atomic system, 
"g" can be engineered to have its desired value using Fes- 
hbach resonance |22j. 

Starting from the action of Eq.Q, the phenomenolog- 
ical GL free energy had been derived by Buzdin et al|23] 
and Combescot et alfMj. In this paper, we will follow 
the notation of Casalbuoni et al[20]. The GL free en- 
ergy functional of a FFLO superconductor in terms of 
the order parameter 0(k) , in Fourier space, is 
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in which a = p[ln(47rr/0o) + Re[-0(l/2 + i5^i/2TTT)]], 
j3 = -(p/167r2T2)Re[?/'(2)(l/2 + iSn/2TTT)] and 7 = 
p/{102ATr'^T^)Re[i;(-^\l/2 + iSn/2nT)]. p{= pI/tt'^vf) is 
the density of states of the electrons at the Fermi surface, 
T is the temperature, ^0 is the superconducting gap at 
zero temperature and zero imbalance, tp{z) and ^'-"•'(z) 
are the digamma function and nth derivative of the 
digamma function, respectively. Here onward we will use 
Eq.([2]) as the starting point of our analysis. We note 
that in Eq.([2]) there is a fourth order gradient term, or 
quartic in momentum, that is quadratic in order parame- 
ter. This term is essential for the occurrence of a Lifshitz 
point. What is vital is that the coefficient of the gradient- 
square, quadratic term and quartic cj) term in Eq.([2| 
changes sign at the same point, namely at /3 = 0. This 
is the LTP which we will take as the hallmark of the GL 
free energy of Eq. ([2| . 

In a general study of tricritical and Lifshitz kind of 
behaviour, Aharony et al[2T] started with the free energy 
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where a denotes an m-dimensional space with m < d 
and subsequently Vq = X^I^i^^^j denotes {d — m) 
dimensional space with = X^iLi™*^*- is an n- 

component order parameter with 0{n) symmetry. In this 
free energy, the momentum dependent terms with cou- 
pling constants Ai and A2 are relevant in addition to w, 



as shown by Aharony et al.|21|. in d < 3 + m/2 dimen- 
sions. This system can have the following multicritical 
points: 

1. Standard tricritical point: t = 0, m = 0, /i 7^ 0; in 
which t is the scaling field related to r, u and /2 
are scaling fields related to u and /i at the critical 
point, t — 0, respectively. 

2. Lifshitz point: i = 0,{( 7^ 0,/2 = 0, i.e. the 
quadratic gradient ^-term and the quartic ^-term 
have simultaneous zeros. 

3. Lifshitz tricritical point: t = 0, u = 0, /i = 0. 

To make it transparent that the GL free energy of FFLO 
state is indeed same as the Lifshitz tricritical free energy, 
we write Eq.([2| in real space 
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The form of the GL free energy in Eq.Q will be needed 
while doing the mean-field analysis. Interestingly, this 
equation also tells us that the FFLO state belong to the 
{d — d), or isotropic, Lifshitz-tricritical class. 

We now analyse Eq.Q to explore the mean- field phase 
diagram of FFLO state. For /3 > 0, the free energy has 
a minimum for a spatially-constant </>, which we call (f>Q, 
and in that case Eq.H reduces to = V[a4>l + I34>q/2 + 
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FIG. 1. Phase diagram of a FFLO superconductor, at mean- 
field level, plotted as a function of a/7 and /?/7. Graphs I and 
II represents the metastable and first-order line, respectively. 
On curve III, the coefficient of the quadratic term in order 
parameter vanishes. 



4 


(a) 




B 


3 

C 

2 




) 




1 






A 1 





.0 0.1 0.2 0.3 


0.4 0.5 0.6 


3.0 
2.5 


(b) 




B 


2.0 








1.0 

o.s 














A 


0.0 


- 3 - 2 - 1 


II 1 2 3 



80 



^40 
t 



20 



(c) 



H(T) 
-21 

-21.5 80 
22 

- 22.5 60 
-23 

24 40 
-28 





22 T 






^1 









2.4 2.6 2.8 3.0 




2.0 



2.5 



3.0 
T(K) 



3.5 



4.0 



FIG. 2. Plot of the specific heat as a function of the scaled 
temperature. Fig. (a) shows as we lower the temperature from 
the system follows specific heat as plotted in curve (1). But at 
the first-order transition temperature, the specific heat jumps 
from point A to B and subsequently follow the mean-field 
specific heat, plotted in curve (2). Fig.(b) shows only the 
jump in specific heat. Fig.(c) is the experimentally obtained 
plot of specific heat as obtained by Lortz et alFTH showing 
the anomalous jump of the specific heat. 



7(/)o/3]. The term proportional to </)q is redundant, and, 
for a < 0, we have the transition from a normal state 
(^0 = 0) to the homogeneous Bardeen-Cooper-Schrieffer 
superconducting state with zero center-of-mass momenta 
Cooper pairs. When (3 turns negative, i.e. when the 
spatially-constant state is no longer the minimum of GL 
free energy, then we envisage a periodic variation of (f) 
written as 0(x) — (p^e^^'^, and with this form of (/>(x), 
the GL free energy in Eq.Q reduces to 
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where V is the volume of the system. Minimising this free 
energy with respect to k and substituting k's minimum 
value in above equation we obtain 



n = v 



5/32 \ 2^ 1370E 
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Since /3 < in Eq.([6]), the minimisation of Eq.Q with 
respect to 00 gives a local minimum of £7 at a finite </>o 
when a = ai — (259/32)/(9367), shown as graph I in 
Fig.([T]), which signals the onset of a metastable state as 
the normal state with 0o = remains the global minima. 
For a lower value of temperature, or a, and for a fixed 
j3 and 7, a first-order transition occurs when the free en- 
ergy at 00 7^ goes below its value at 0o = 0. This 
occurs at a = a2 = (27/32)/(1047), shown as graph II in 
Fig.g. At even lower temperature, the coefficient of the 
quadratic (/>o-term in Eq.([6]) goes to zero. The tempera- 
ture for this line is given by a = = (5/32)/(247) and 
this has been shown as graph III in Fig. ([!]). 

With this phase diagram available, we are now 
ready to find the behaviour of specific heat near the 
first-order transition. Above the first-order transition 
temperature (02), the mean- field specific heat is iden- 
tically zero. To calculate its value above the first- 
order line, we consider gaussian fluctuations around 
the mean field value and calculate the partition func- 
tion Z = J l'0exp[— r^g/fcsTc"] in which Q.g — 
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s the gaussian free en- 
ergy. We have considered /3 < as we are below the 
metastable line. We minimise this free energy with re- 
spect to k and write k = kp + 1 where fc^ = 5|/3|/87. 
With this partition function, we calculate the thermo- 
dynamic free energy per unit volume which \s F = 
—ksTc^nZ and, taking the second derivative of F with 
respect to a, the specific heat per unit volume. The form 
of the gaussian specific heat, which is plotted as curve 
(1) in Fig. (2a), comes out to be 
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Here we have assumed that the fluctuations around kp 
are three dimensional. That is why our exponent is differ- 
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ent from that of Konschelle et al[5S] who took the fluc- 
tuations in the radial direction alone. We further ob- 
serve that inclusion of the higher order terms will lead 
to a substantial increase in the rate of divergence of Cg 
as shown by Aharony et al[lT] who evaluated the first- 
order correction to the specific heat using renormaliza- 
tion group techniques. However, for temperatures above 
the first-order temperature, which is greater than the 
temperature at which Cg diverges, the gaussian specific 
heat is an accurate representation of the fluctuation spe- 
cific heat. To find the specific heat of the system below 
the first-order temperature, one needs to evaluate the 
mean-field specific heat. This can be calculated by tak- 
ing the double derivative of the mean-field GL free energy 
of Eq.Q with respect to a, the temperature parameter. 
The mean-field specific heat per unit volume, plotted as 
curve (2) in Fig. (2a), is 



(8) 



With these forms available, let us now discuss the 
experimental picture as one measures the specific heat. 
When one decreases the temperature of the imbalanced 
fermionic system, the relevant specific heat that one mea- 
sures is the gaussian specific heat. But then the system 
confronts the first-order line and the specific heat jumps 
to its mean-field value. This, we claim, is the source of 



anomalous jump in the specific heat as observed by Lortz 
et al in their experiment. Considering the experimental 
data for specific heat jump at an magnetic field, applied 
parallel to the superconducting layers, of 22 T, we fix 
(a/3)^/^ = 0.025. The jump of the specific heat is plot- 
ted in Fig. (2b) and we have also shown the experimental 
specific heat diagram from Lortz et al in Fig. (2c). Thus 
we have shown that the static analysis of the GL free 
energy provides the reason for the jump in specific heat. 
Now one can legitimately ask how to pinpoint the region 
in the phase diagram where one will find the jump in 
specific heat in an experiment? 

To provide an answer to this question, we have studied 
the dynamics of this system near LTP using TDGL for- 
malism. This study has also been encouraged by a very 
recent study by Sodemann et al|26| who showed that 
finite-momentum superfluid state occupies a large area 
in the phase diagram of an imbalanced fermionic system 
just following a quench from the normal to the supercon- 
ducting state than in the equilibrium phase diagram of 
the same system. In our study, we look for the growth 
of FFLO phase following a quench. The TDGL for this 
system can be written as dcf)k/dt — ~r{Sil/54>^k), in 
Fourier space where Q is the GL free energy in Eq.([2|. 
The quantity of experimental relevance is the dynami- 
cal correlation function or the structure factor- S'(k, t) = 
(^k(i)</'-k(0)- Using the TDGL we get 
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To make progress, we perform the Hartree approxima- 
tion, in which the correlation functions factor, to the 
quartic and sextic correlations. Implementing this on 
Eq.Q, we arrive at 

55k 
dt 
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From this equation it follows that at very short time, 
when linear approximation holds, all modes decay for /3 < 
if we are in the region II of the phase diagram(Fig.([T])). 
To look for the long-time dynamics we work, as before, 
for /3 < and expand the wave number around k^, i.e. 



write k 



1. Restricting to terms quadratic in T 



we can rewrite Eq.(10l as 
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This allows us to write the solution of Eq.( 11 ) as 



Sk{t) = Acxp 
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where g{t) = £dt'f{t'). If one now demand, as one 
does in this kind of self-consistent calculation [27l |28] . 
that J2p g06S to some constant value as t — )■ oo, thus 
obtaining the ordered state; then for t — > cx) one can ne- 
glect the sub-leading J2pP^'^p term, and, the final result 
for the structure factor comes out to be 
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This result shows that at long time after the quench 
5'k(i) — > for all wave numbers expect k = kc. So de- 
termination of the parameter fc^ actually fixes the ratio 
/3/7 of the phenomenological model of Eq.([2]). 

In conclusion our analysis provides a theoretical ex- 
planation for the experimental observation of the FFLO 
state in layered organic superconductors. In particular, 
we have studied the FFLO phase diagram using the phe- 
nomenological GL free energy and found that the phase 
diagram contains a LTP. We have calculated the specific 
heat and found that it jumps near the first-order line in 
the phase diagram as we lower the temperature, as also 
found in an experiment by Lortz et al|14j. By doing a 
dynamical analysis of the GL free energy, we suggest a 
way to explore the portion of the phase diagram where 
one should look for the jump. A possible extension of 
our work can be performing the dynamical analysis be- 
yond the Hartree approximation and exploring the elu- 
sive FFLO state further. 
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